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Abstract 

We prove the following conjecture of Zeilberger. Denoting by C n the Catalan 
number, define inductively A n by (-l)"" 1 ^ = C n + Y%Zl(-l) j C n -j and 

a n = 2A n /C n . Then a n (hence A n ) is a positive integer. 



1 Introduction 

Let z be an indeterminate. For any nonnegative integer r, we denote by C r (z) the 
Narayana polynomial denned by Cq(z) = 1 and 



C r (z) = J2N(r, k)z k ~ x 

k=l 

with the Narayana numbers N(r, k) given by 

N(r, k) 



r \k — 1 J \k 

There is a rich combinatorial litterature on this subject. Here we shall only refer to [T4~| 
T8] , fTf\ Exercice 36] and references therein. 

We have C r ([T]) = C r , the ordinary Catalan number, since 



v 7 k=i 



Moreover [1] the Narayana polynomial C r (z) can be expressed in terms of Catalan numbers 



as 



C r (z) = J2z m (z + 1 ) r - 2m ~ 1 ( r 2 ^)Cm. (1) 
m>0 ^ 7 
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The polynomial 

(z + 1) C r (z) - C r+1 (z) = (1 - r)z r - 1 + ... 

can be expressed, in a unique way, in terms of the monic polynomials z m C r -2m+i( z ) with 
degree r — m. As a consequence, we may define the real numbers A m (r) by 



(z + 1) c r (z) - c r+1 ( Z ) = j2(- z r (L _\) 

m>l ^ ' 



)• (2) 



Note that the real numbers {A m (r),m > 1} depend on r and that we have A\{r) = 1. 
By relation ([!]), the left-hand side may be written as 



fc>0 

Similarly the right-hand side becomes 



£ (-ir*-n* + i) r - 2m - 2n ( 2 r m _\) ( r 2n 2m W r) °»- 



m>l,n>0 

Comparing coefficients of z k (z + l) r ~ 2fc on both sides yields 

\ 7 m+n=k v 7 v 7 

with r > 2&. This is easily transformed to 



Cn = i2(-iy- 1 (^ j _ 1 ^Mr)c^. 



In other words, A n := A n (r) is an integer independent of r, given by the recurrence 
formula 

n— 1 

. r>n — | \ 

(3) 

This relation can be chosen as a definition, equivalent to 
In Section [2] we prove 



Theorem 1. T/ie integers {A n ,n > 2} are positive and increasing. 
The values of A„ for 1 < n < 14 are given by 

1, 1, 5, 56, 1092, 32670, 1387815, 79389310, 5882844968, 548129834616, 62720089624920, 
8646340208462880, 1413380381699497200, 270316008395632253340. 
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At the time the positivity of A n was only conjectured, we asked Doron Zeilberger for 
an advice. He suggested to consider also the numbers a n = 2A n /C n which, by an easy 
transformation of (jSJ), are inductively defined by 



In Section [3] we prove the following conjecture of Zeilberger. 

Theorem 2. The numbers {a n ,n > 2} are increasing positive integers. They are odd if 
and only if n = 2 k — 2 for some k > 2. 

The values of a n for 1 < n < 16 are given by 

2, 1, 2, 8, 52, 495, 6470, 111034, 2419928, 65269092, 2133844440, 83133090480, 
3805035352536, 202147745618247, 12336516593999598, 857054350280418290. 

Both sequences seem to be new. 

Of course since 2A n = a n C n , Theorem [1] is an obvious consequence of Theorem [2j 
However we begin by a direct proof of Theorem [TJ which has its own interest. This proof 
describes the algebraic framework of our method (the theory of symmetric functions), and 
it yields a generating function for the A^s. 

The An's are also worth a separate study in view of their connection with probability 
theory, which was recently noticed. Actually Novak [15J observed, as an empirical evi- 
dence, that the integers (— l) n ~ l A n are precisely the (classical) cumulants of a standard 
semicircular random variable. Independently Josuat- Verges [9] defined a g-semicircular 
law, which specializes to the standard semicircular law at q — 0. He described combinato- 
rial properties of its (classical) cumulants, which are polynomials in q having (— l) n ~ l A n 
as their constant terms. It follows from his work that A n may be obtained by a weighted 
enumeration of connected matchings (equivalently, fixed-point free involutions on 1, 2n 
such that no proper interval is stable). Finding a similar bijective proof for a n would be 
very interesting. 

2 Properties of A n 

This section is devoted to a proof of Theorem [IJ We shall use two remarks of Krattenthaler 
and Lascoux. Define 




n — j + 1 
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Krattenthaler observed that the definition ([3]) is equivalent with 



A(z)C(z) = 2^C(z). (4) 



Therefore we have only to show that, if we write 

d 
dz 



C'(z)/C(z) = £log(C(z)) = £ 



r ~m— 1 



m>l 

the coefficient c m has sign (— l) m_1 . 

Let § denote the ring of symmetric functions [321 Section 1.2]. Consider the classical 
bases of complete functions h n and power sums p n . Denote 

n>0 n>l 

their generating functions. It is well known [T31 p. 23] that P(z) = H'(z)/H(z). 

Since the complete symmetric functions h n are algebraically independent, they may 
be specialized in any way. More precisely, for any sequence of numbers {c n ,n > 0} with 
Co = 1, there is a homomorphism from § into the ring of real numbers, taking h n into 
c n . Under the extension of this ring homomorphism to formal power series, the image of 
H(z) is the generating function of the c n 's. By abuse of notation, we write h n = c n and 

Therefore it is sufficient to prove the following statement. 
Lemma 1. When the complete symmetric functions are specialized to 

1 



n\(n + 1)! 

i.e. when H(z) = C(z), the coefficients of P(—z) are all positive. 

Lascoux observed an empirical evidence for the following more general result (which 
gives Lemma [1] when x — 2). 

Lemma 2. Let x be a positive real number and (x)k = Yli=i( x ^ ~ -0 denote the classical 
rising factorial. When the complete symmetric functions are specialized to 



n\(x) n 

the coefficients of P(—z) are all positive. 
Proof. With this specialization we have 

1 z n 

n>0 
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the so called confluent hypergeometric limit function. Since 

— Fi{x; z) = -oF^x + l;z), 
az x 

we have 

oFAx + l-z) 



P z) 



xqF^x; z) 

The classical Gauss's continued fraction [201 P- 347] gives the following expression for the 
right-hand side 

Mx + l-z) 1 



XqF^x; z) 



{x + 1) + 



(x + 2) + 



(x + 3) + '•• 



This continued fraction may be written as a Taylor series by iterating the usual binomial 
formula 



-2 



fe 2 >0 

where /f 2 is itself of the form (ci(l + f 2 z))~ k2 . This classical method ( [T2~t Exercise 
4.2], [16J) yields 

1 i-r / -z \ k ' fa + k l+1 - 1 

(kiMM,-) *>i v v t+i 

where the fcj's are nonnegative integers. This proves the lemma. □ 

In the previous formula, the contributions with fc, = and fci+i > are necessarily 
zero. Hence for n > 2 we have 



e n 



with C n the set of sequences k = (ki, k 2 , . . . , k{) of I < n — 1 (strictly) positive integers 
summing to n — 1. There are 2 n_2 such sequences (compositions of n — 1). For instance 

1 _ -1 _ 1 1 

^ x' ^ 2 x 2 (x + l)' ^ 3 x 2 (x + l) 2 (x + 2) x 3 (x + l) 2 

To each k = (ki,k 2 , . . . ,k t ) e C n let us associate two elements of C n+ i defined by 
«4 = (k%, k 2 , . . . , ki + 1) and k 2 = k 2 , . . . , h, 1). We have 

C„+i = [J {ki, k 2 }. 



Denoting g K ,g Kl ,g K2 the respective contributions to (j5J) of k, ki, k 2 , we have easily 

-(x + n - l)(x + n)g Ki > g K (i = 1,2), 

the equality occuring for k = (1, 1, . . . , 1) and i — 2. Summing all contributions to (JSJ), 
we obtain 

-(x + n - l)(x + n)p n+1 > 2p n . 

Since (jlj) implies 

A n = (-l) n - 1 2(2n-iy. Pn \ x=2 , (6) 

substituting x = 2, the above estimate yields A n+ i > A n , which achieves the proof of 
Theorem [TJ □ 
Incidentally we have shown that the generating function of the A n 's is given by 

A 9 
1 ' (2m -1)! z 

m>l v ' o _j_ 



3 + 



z 

4+ 



5+ '•• 

3 Properties of a n 

This section is devoted to the proof of Theorem [2j The following lemma plays a crucial 
role. 

Lemma 3. When the complete symmetric functions are specialized to h n = l/((x) n n!) 
with x > 0, we have 

zP(z) 2 + z^- P{z) + xP{z) = 1. 
Equivalently for any integer n > 1, we have 

n 

(n + X)p n+ i = -^PrPn-r+l- (7) 
r=l 

Proof. In view of 

P(z) - o F i(x + i;z) 



x Fi(x; z) 
we have 

d 

dz y ' x 2 



P(z) 2 +^-P(z) = -—^ -( F 1 (x+l;z) 2 +^— F 1 (x; *) Fi(x+2; *)- Fi(x+1; zf 

~qPi{x;z) z \ x + 1 



hence 



d n/ N i __ nf ^ z F 1 (x + 2;z) i o-Pi + 1; -z) 



zP(z) J + z-j- P(z) + zP(z) = — " . \ J + 



dz x(x + 1) Fi (x; z) Fi(x;z) 
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We have 



since 



3 Fi(x; z) - F 1 {x + l]z) = ———oF^x + 2; z) 

x[x + 1) 



~fe / 1 I \ -V -vfc— 1 



fc! \(a:)fc (x + l) k J x(x + 1) (A; - 1)! (x + 2) fc _j 

Substituting x = 2, we obtain a second (inductive) proof of Theorem HJ 
Theorem 3. The integers {A n ,n > 2} are positive, increasing and given by 



□ 



. vA n - r + 1 (2n + 1 \ , , 



n + 2 V2r-1 

r=l v 



T/ie numbers {a n ,n > 2} are positive, increasing and given by 

n + l\ fn + 1 



lA 1 /n + 1 
2^n + l\r+l 

r=l x 



Proof. We substitute ([6]) into (J7J). Then we substitute A n = a n C n /2. In both cases we 
proceed by induction. □ 



The positive integer 



n + 1\ / n + 1 \ / n \ / n\ ( n \ ( n 



n + 1 \r + 1/ \r — 1/ \r — 1/ \r/ \r — 2/ \r + 1 

is known [7J [8] to be the number of walks of n unit steps on the square lattice (i.e. either 
up, down, right or left) starting from (0,0), finishing at (2r — n — 1,1) and remaining in 
the upper half-plane y > 0. One has c n r = c n n _ r+ i, which corresponds to symmetry with 
respect to the y-axis. 

We need some arithmetic results about c n r . Given a rational number r, there exist 
unique integers a,b,m with a and b odd such that r = 2 m a/b. The integer m is called 
the 2-adic valuation of r and denoted ^(r). A classical theorem of Kummer (appeared 
in 1852, see [TH] for an inductive proof) states that the 2-adic valuation of ( a ~^ b ) is a 
nonnegative integer, equal to the number of "carries" necessary to add a and b in base 2. 

Lemma 4. (i) If n is even, c n>r is even. 

(ii) If n is odd and r even, c n , r is a multiple of A. 

Proof, (i) Since n + 1 is odd, we have 

Mc n ,r) = 1 + VI ( ( n + }))+V2( (" + }\ ) > 1. 



r + 1/ / \\r — 1 



7 



(ii) Since r + 1 is odd, we have 



Z*2(C„, r ) = 1 + Z/ 2 I [})+V2 " 



r I I V V r — 1 



Since r — 1 and n — r + 2 are odd, at least one carry is needed to add them in base 2. 
Thusz/^ 1 ))^!. □ 



We also consider 

1 ( 2p N 2 
C2p ~ 1 ' p " p U + 1 



Lemma 5. (%) 7/p zs even, c 2f ,_i iP zs a multiple of 8. 

(ii) If p is odd and p ^ 2 k — 1 for any k > 1, c 2p _i )P is a multiple of 4. 
(m,) 7/p = 2 fc — 1 /or some fe > 1, c 2p _i iP is odd. 



Proof, (i) We have 



C2p-l,p 



Ap (2p-V 2 



(p + l) 2 V p 



Since p + 1 is odd, we have 

^2(c2 P -i,p) = 3 + 2z/ 2 (J^ P p J J - 3 - 
Since p is odd, we have 

We may write the binary representations of p + 1 and p — 1 as 

p + 1 = c k+m c k+m -i . . . c fc+2 10 ... 00, 
p - 1 = c fe+m c fc+m _i . . . c fc+2 01 ... 10. 

At least one carry is needed to add them, except if Ck+ m = Ck+m-i = ■ ■ ■ = c/ c+2 = 0, i.e. 
p + 1 = 2 fc , which is excluded. Therefore ^((^J) > 1- 

(Hi) When p + 1 = 2 k , the binary representations of p + 1 and p — 1 are 10 ... 00 and 
01 . . . 10. Therefore ^((^J ) = 0. □ 

Our proof of Theorem [2] also needs the following result. 

Theorem 4. T/ie numbers {A n ,n > 3} and {C n ,n > 3} /iai>e t/ie same parity. They are 
odd if and only if n = 2 k — 1 /or some k > 2. 
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Proof. It is known [5j [TT] that the Catalan numbers {C n ,n > 3} are odd if and only if 
n = 2 k — 1 for some k > 2. Therefore we have only to prove the first statement. 

We proceed by induction on n. We use ([2]) to express (— l) n A n + C n , as a signed sum 
of the terms 

(2n - 1\ 
% = V2J-V j 
for 1 < j < n — 1. Two cases must be considered. 

Case 1: n = 2 m for some m. Then yi and y n _i are obviously odd. On the other hand, 
the contributions {yj, 2 < j < n — 2} are even, because the conditions n = 2 m , j = 2 k — 1 
and n — j = 2 l — 1 cannot be simultaneously satisfied if j ^ 1 or n — j ^ 1. 

Case 2: n^2 m for any m. Then all contributions {yj, 1 < j ; < n — 1} are even. Indeed 
if and Cn-j are odd, we have j = 2 k — 1 and n — j = 2 l — 1 with k, I > 2, since j and 
n — j cannot be equal to 1. But then (^Zj.) is even because we have 



2n - 1\ _ /a + 6' 
2.7 -V V a . 

with a = 2 fe+1 — 3 and b = 2 l+1 — 2. Then the binary representations of a and b are 
1 . . . 1101 (k + 1 digits) and 1 . . . 1110 (/ + 1 digits). Since k,l > 2, at least one carry is 
needed to add them. Summing all contributions, A n and C n have the same parity. □ 

Proof of Theorem \M We proceed by induction on n. Our inductive hypothesis is 
twofold. We assume that for any m < n, the positive number a m is an integer, and 
that a m is odd if and only if m = 2 k — 2 for some k > 2. 

If n is even, say n = 2p, a n+ i is an integer because by Theorem 3 we have 



1 P 

a-2 P +i = - y c 2p 



v 

r Qj r Qj2p—r+li 



2 

r=l 

and we apply Lemma |4] (i). Moreover a 2p +i is even, since a r or a 2p _ r +i is even for each r. 
If n is odd, say n = 2p + 1, we have 

l ^ > l 2 /Q\ 

a 2p+2 — — / J C2p+l,r a r a 2p-r+2 + ^C2p+l,p+l a p +i- (8) 
r=l 

Each term of the form C2 P +i. r or a r a2 P ~ r +2 is an even integer. This results from Lemma H] 
(ii) if r is even, and from the fact that a r and a>2p- r +2 are even, if r is odd. We are left to 
consider the term 

1 

W p — ^ C 2p+l,p+l a p+l- 

This is an integer whenever a p+1 is even. If a p+1 is odd, by our induction hypothesis, one 
has p+ 1 = 2 k — 2 for some k > 2, and we apply Lemma [5] (i). Summing all contributions, 
&2p+2 is an integer. 
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It remains to show that a2 P +2 is odd if and only if 2p + 2 = 2 k — 2, i.e. p + 1 = 2 k ~ 1 — 1 , 
for some k. We have seen that all terms on the right hand side of (|S]) are even, except 
possibly for w p . 

If p+1 is even, w p is even because of Lemma [5] (i). If p+1 is odd and p + 1 ^ 2 k — 1 for 
any k, w p is even by Lemma [5] (ii) and because a p+ i is even. Finally the only possibility 
for w p to be odd occurs when p + 1 = 2 k — 1 for some k. Then in view of Theorem HI 
A p+ i and C p+ i are odd. Hence Op+i/4 = (A p+ i/C p+ i) 2 is also odd. Applying Lemma [5] 
(iii), w p is odd. □ 



4 Extension to type B 

Generalized Narayana numbers have been introduced in [HI Section 5.2] in the context of 
the non-crossing partition lattice for the reflection group associated with a root system. 
Ordinary Narayana polynomials correspond to a root system of type A. For a root system 
of type B, generalized Narayana polynomials are defined [6, Example 5.8] by Wq(z) = 1 
and 

fc=0 v 7 

For their combinatorial study we refer to [2| |3] and references therein. 
We have W 7 ,.© = W r , the central binomial coefficient, since 

v 7 fc=0 v 

Moreover [3l equation (2.1)] the Narayana polynomial W r (z) may be expressed in terms 
of central binomial coefficients as 



2 



(9) 



W r (z) = Y,Z m ( Z + 1 Y~ 2m L r m )W m . 
m>0 ^ 7 

Exactly as in Section [2] the polynomial 

(z + 1) W r (z) - W r+1 (z) = -2rz r + ... 

may be expressed, in a unique way, in terms of the monic polynomials z m W r -2m+i{z) with 
degree r — m + 1. As a consequence, we may define the real numbers B m (r) by 



2m+l{z). (10) 



(z + 1) W r (z) - W r+1 {z) = J2(- Z ) m ( 2m r _ ! WW W r . 

m>l ^ 7 

Note that the real numbers {B m (r),m > 1} depend on r and that we have Bi(r) = 2. 

Using relation (Q, and paralleling the computation in Section [H this definition is easily 
shown to be equivalent to 



B 3 {r)W n . r 
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Therefore B n := B n (r) is an integer independent of r, given by the recurrence formula 

71-1 

~3 1 



n-1 _ 1 \ 

•l)- 1 ^ = iy n + [J. _ i )b j w n . 



id 



which can be chosen as a definition, equivalent to ( flOl) . 
As in Section |2} if we define 

n>0 v 7 n>0 v 7 

B( z ) = 5m ~— i 



m>l 

the definition ( fill) is equivalent with 



(2m- 1)! 



B(z)W(z) = 2^-W{z). 

Obviously W(z) = o-fi(l;#)- Therefore we are in the situation of Lemma |2] with x = 1 
and we have 

B n = (-l) n - 1 2(2n-l)!p n U =1 . (12) 

As a consequence, the integers {B n ,n > 1} are positive and increasing. Moreover 
their generating function is given by 

R 2 
1 ' (2m -1)! z 

m>l v ' 1 j 



Z 

3 + 



4+ ••• 

By an easy transformation of ffTTj) . it is clear that the positive numbers b n = B n /W n 
are inductively defined by 

71-1 



-^-i + d-*(;::)G)*. 



which proves that 6 n is an integer. This sequence of positive integers is not new. It was 
studied by Carlitz [TJ equation fl3])] in connection with Bessel functions. 
Substituting (JT2]) into Lemma [3j we have also 



r=l x ' 

bn+1 = ( ) ( i W&n-r+l- 



r-'-'n— r+1 ? 
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The last equation was already mentioned in pQ. It is very similar to Theorem [31 Actually 
the positive integer 

/ n\ ( n 

1, 

is known [7J [S] to be the number of walks of n unit steps on the square lattice (i.e. either 
up, down, right or left) starting from (0,0) and finishing at (2r — n — 1, 1). 

Remark: The Narayana polynomial of type D is given [61 Figure 5.12] by W r (z) — 
rzC r -\(z). There is an analog of fl2]) and (10), but it does not seem to involve an in- 
teger sequence. 



5 Open problems 

As indicated at the end of the introduction, we are in lack of a combinatorial interpretation 
for a n , possibly in the framework of probability theory. Here are two open questions which 
might be also investigated. 



5.1 Schur functions 

It is known (T3J p. 23] that the elementary symmetric functions e n satisfy 

n 

Tie n ^ ^ ( 1) Pr&n—r- 
r=l 

Lemma [2] states that, when h n is specialized to l/((x) n n\) with x > 0, (— l) n_1 pn is 
positive. By an easy induction this implies the positivity of e n . 

Since h n = s n and e n = S]n are the row and column Schur functions, we are led to the 
following conjecture (checked for any partition A with weight < 12). 

Conjecture. When the complete symmetric functions are specialized to h n = l/((x) n n\) 
with x > 0, any Schur function s\ takes positive values. 

A nice expression of s\ in terms of x may also exist. For elementary symmetric func- 
tions e n , such a formula can be obtained as a consequence of their generating function [2T] 



E(Z) = Y,CnZ n = —, r = 



u 2 z 

1 - u 2 z H 

U4Z 

1 - u 3 z H 

1 — u^z + ' • ■ 



with u n = l/(n(x + n — 1)). For another way of writing E( continuous fraction 

involving ratios of Schur functions, we refer to [121 Sections 4.3 and 5.2]. 
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5.2 An identity 

By a classical result [TO], p. 335] we have 

2 / 2r 



Cr+l — 2 C r 

r \r — 2 



For z — 1 the definition (T2J) yields 

2 C n — C n+ i = ^ ^ ( — l) m ( j^ro Cn-2m+lj 



m>l 

hence 



m>l v 7 



2 / 2n 

•2m+l — I 

n\n — 2 



This formula can be used in two ways to define A n inductively. Actually if we write it for 
n = 2r and n = 2r + 1, we obtain 



4r \ ^-A , Sm / 2r — 1 



m=l 

r-1 



'm ( ~ / 2r— 2m+l; 



T(i r >E<- i K 2m 2r - 1 )^ 

\ / m=1 \ / 



2r + 

Both definitions of A r must be consistent, which yields 



2r-2m+2- 



m=l \ \ ✓ \ / 



4r \ 2 /4r + 2 



2r-2/ 2r + lV2r-l 



Dividing both sides by C 2r we obtain 



r-1 



E, l\ m A m r — m y-r (2r — i) 2 2r(r — 1] 

V A.) (9m - 1\\ Or - Om X 3 11 



4/ (2m- 1)! 2r-2m + 3 A_l 4 r _2i- 1 (2r + l)(2r + 3) 

a very strange identity for which we are in lack of an interpretation. 

A similar identity is obtained for the B n 's by using relation ffTUj) for z — 1 and paral- 
leling the above computation. We get 

§ ( - irB ™ (( 2m 2 - 1) w ^ « - 4r Gm :\) ) 

= (4r + 2)W 2r - 8riy 2r _i - W 2r+ i. 
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Dividing both sides by W 2r this yields 

^/ l\ m B m (r - m) 2 2 pj- 1 (2r - i) 2 r— 1 

^ V 4/ (2m - 1)! (2r - 2m + l) 2 -1-1 4r - 2z - 1 ~ ~2(2r + 1)' 

m=l \ / \ / \ / 

In view of ([6]) and (1121) it would be interesting to investigate whether such an identity 
is a general property of the p n 's. 
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